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ABSTRACT 

We report the results of the first transit timing variation (TTV) analysis of the 
very hot Jupiter OGLE-TR-132b, using ten transits collected over a seven-year 
period. Our analysis combines three previously published transit light curves 
with seven new transits, which were observed between February 2008 and May 
2009 with the new MagIC-e2V instrument on the Magellan Telescopes in Chile. 
We provide a revised planetary radius of = 1.23 it 0.07 Rj, which is slightly 
larger than, but consistent within the errors, of the previously published results. 
Analysis of the planet-to-star radius ratio, orbital separation, inclination and 
transit duration reveals no apparent variation in any of those parameters during 
the time span observed. We also find no sign of transit timing variations larger 
than —108 ±49 s, with most residuals very close to zero. This allows us to place 
an upper limit of 5-10 for a coplanar, low-eccentricity perturber in either the 
2:1 or 3:2 mean-motion resonance with OGLE-TR-132b. We similarly find that 
the data are entirely consistent with a constant orbital period and there is no 
evidence for orbital decay within the limits of precision of our data. 
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Introduction 



Precise timing of exoplanet transits over many years is a powerful technique to learn 
more about a planetary system. Transit timing variations, or TTVs, can arise from many 
potential scenarios, such as orbital precession, orbital decay, and perturbations by additional 
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Recently, iHolman et al.l (120101 ) announced the first incontrovertible evidence of transit 
timing variations in the Kepler-9 system, where a pair of Saturn-mass planets are trapped 
in a 2:1 mean- motion resonance, producing strong, complementary deviations of over an 
hour in the midtimes of both transits. Other observational evidence is less clear cut, with 



recent claims of timing v ariations in WASP-3b ( iMaciejewski et al.ll2010bl ) and WASP-lOb 



(Adams et al. 


2010a; 


Diaz et al. 


2008) 



there have been claims of a tentative detection of variations in the transit duration of G J436b, 



increasing by roughly 3 minutes per year ( ICoughlin et al.ll2008l ). con sistent with the < 8 
limit for a second planet in this system placed by transit timing (IBean k SeifahrtI l2008l ). 
Finally, previous work by our group has found very tentative evidence that the orbital 
period of 0GLE-TR-113b is decreasing at a r ate of —60 =1= 15 ms yr~^, though additional 
observations are needed to confirm this result ( lAdams et al.ll2010b( ). 



Though perhaps less exciting, the non-detection of transit timing variations is both im- 
portant and expected. It has been suggested t hat systems with h ot Jupiters are statistically 
less likely to be part of multi-planet systems ( IWright et al.ll2009l ). and transit timing limits 
can place tighter mass limits on potential nearby companions than other detection methods 
to date, with upper limits of a few Earth masses possible for companions in strong mean- 
motion resonances. With observations spanning many years, upper limits can also be placed 
on the rate of orbital decay of the known hot Jupiter, providing valuable direct observational 
constraints on theoretical models of tidal dissipation in stars. Currently tidal dissipation in 
stars is only estimated indirect ly, through arguments about circularization timescales of 
binary stars in clusters (see e.g. iMeibom k Mathieul 120051 ; lOgilvie k Linl 120071 ). 



In this work we present the case of the hot Jupiter OGLE-TR-132b, with seven new 
light curves observed and no timing variations larger than — 108 ± 49 s detected during 
observations spanning seven years. In Section |2] we describe the observations, data analysis 
and light curve fitting. Section 12] describes the modeling and improvements to the system 
parameters, and places limits on parameter variability between transits. Section H] discusses 
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the timing constraints placed on the planetary system. We conclude in Section [51 



2. Observations and photometry 



OGLE-TR-132b was first announced as a planet candida te by lUdalski et al.l (120031 ) 
and confirmed to be of planetary mass by iBouchy et al.l (120041 ) . The star, OGLE-TR-132, 
is one of the faintest transit-hosting stars (1=15.7; RA(J2000)=10:50:34.72, Dec(J2000)=- 
61:57:25.9). It is somewhat larger and hotter than the sun, with = 1.26 ±0.03 Mq, i?^, = 
1.34±0.08 Rcr,. T,tf = 6210±59ir, and relatively high metalhcity, [Fe/H] = +0.3 7±0.07 dex 



iBouchv et al.1 120041 ). No reliable estimate of the age of OGLE-TR-132 exis ts (iMelo et al. 



2006h . Tw o high-quality follo w-up light curves in i?-band were obtained by iMoutou et al.l 



( I2OO4J ) and iGillon et al.l ( 120071 ) , and remain until this work the only published photometry of 



this p lanet aside from the discovery data. Comb i ning t hose transits with the iBouchy et al. 



( I2OO4I ) radial velocity observations, iGillon et al.l (120071 ) derived a radius for this planet of 



Rp — 
of R, 



1.18 ±0.07 Rj, which is slightly larger, but consistent within the er rors, than th e value 
= 1.13 ± 0.08 Rj derived by Moutou et el. (2004). More recently, ISouthworthI J2OO8I ) 



reanalyzed the available data and derived a new radius for the planet, Rp = 1.25 ± 0.16 Rj, 
and identified this system as one that needs additional photometry to improve the light-curve 
parameters. 



Since 2006 we have conducted a campaign to search for transit timing variations ([Adams 



2OIOI ). using the same instrument and observing strategy to minimize systematic differences 
between d ata sets. Results for OGLE-TR-lllb and 0GLE-TR-113b have been previously 
published ( lAdams et al.ll2010al 1bl). OGLE-TR-132 is our faintest target to date, but because 
of the short period of the planet (P = 1.69 days) there are many transit opportunities. We 
observed seven full transits of OGLE-TR-132b between April 2008 and May 2009, using 
the dual-CCD instrument MagIC on the Magellan-Baade 6.5m telescope. All transits were 
observed in the Sloan i' filter using MagIC-e2v; an eighth, partial transit was attempted 
with MaglC-SiTe in February 2008, but a faulty shutter rendered those data unusable. In 
the remainder of the paper we refer to each transit by the UT date when it was observed, 
following the notation YYYYMMDD. 

The MagIC-e2v CCD saw first light on Magellan in April 2008. The camera has low 
readout noise (about 6 e- per pixel), a small field-of-view (38" x 38") and high- resolution 
pixels (0."037 per pixel), which both minimizes blends and produces images of the stars 
spread over many pixels (FWHM=6-20 pixels, depending on seeing and binning). This 
latter feature reduces differential pixel response effects and increases the total number of 
photons that can be collected per frame, without requiring defocusing. The e2v gain was 2.4 
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e-/ADU in 2008 and 0.5 e-/ADU in 2009 (after engineering changes). The main advantage 
of the e2v over conventional CCDs is its frame-transfer capabihty: the readout time per 
frame is only 5s in standard readout mode and a few milliseconds in frame transfer mode, 
surpassing the 23s readout time of the MaglC-SITe chip. 

Given that the main focus of our observations is accurate timing, we took special care 
to correctly record the time in the image headers. The 2008 MagIC-e2v times came from 
the network server, which we verified to be synchronized with the observatory's GPS clocks 
at the beginning of each night. The 2009 MagIC-e2v times came from an embedded PC104 
computer, which receives unlabeled GPS pulses every second, and is synchronized to the 
observatory's GPS receiver every night. The intrinsic error for all header times is thus much 
less than a second. 

The observational conditions for each transit are summarized in Table [TJ Because 
MagIC-e2v is a new camera, we have observed with several different configurations in order 
to test performance. All three transits in 2008 were observed in the same configuration: 5s 
readout-mode and 1x1 binning. Three transits in 2009 used 2x2 binning, while one transit 
(20090424) was binned 1x1; all 2009 transits used the frame-transfer mode with millisecond 
readout times. Exposure times ranged from 10-120 s depending on the binning and the 
weather conditions, and were adjusted with a goal of having at least 10^ counts on the target 
on each frame. A few images were discarded because of low counts or telescope tracking 
problems during the observations, as noted in Table [H For calibration data, in 2008 we used 
twilight flats (5 per night), while in 2009 we switched to dome fiats because we could get 
many more of them (typically 100 per night), thus reducing the potential noise introduced 
during the calibration stage. Although bias frames were taken each night, we did not use 
them, preferring instead to subtract the bias level recorded in the overscan region on each 
data frame. 

The data for each night were overscan-corrected and fiatfielded using standard IRAF 
routine^. Simple aperture photometry, using the IRAF apphot package, provided optimal 
results, since both the target and several comparison stars appear well isolated in all the 
frames. We ran a wide grid of apertures and sky annuli to find the values which provided 
the optimal photometry, i.e. the values which minimized the out-of-transit photometric 
dispersion of each light curve. The comparison stars were iteratively selected to be 1) similar 
in brightness to the target, 2) photometrically stable, and 3) un-blended. Of the 10-20 



^IRAF is distributed by the National Optical Astronomy Observatory, operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foun- 
dation 
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stars examined for each transit, we selected between 2-8 as final photometric comparisons, 
depending on the night. The best aperture radii ranged between 9-26 pixels, while the best 
sky annuli had inner radii between 20-50 pixels and 10-30 pixel widths. The best photometry 
parameters are compiled in Table [2J 

The final step in the photometric analysis is the identification of potential systematic 
trends against all available parameters. We examined the out-of-transit baselines of each light 
curve by fitting linear trends with respect to airmass, seeing, telescope azimuth, (x,y) CCD 
pixel location, and time. The 20080419 transit light curve showed no significant trends, while 
the other light curves showed trends against telescope azimuth (20080424), seeing (20090424), 
and airmass (20080511, 20090207, 2009031, 20090511). With each trend removed, there is a 
possibility of introducing error into the light curve, though we note that three of our light 
curves are within 20% of the photon noise limit. 

We now briefiy describe the previously published light curves that we used in our analy- 
sis. All ten transit light curves are shown in Figured! where we have averaged the data in 120 
s bins for representational clarity. The full unbinned light curves are available electronically, 
and are excerpted in Table El 



2.1. Literature light curves 



Two i?-band light curves were reported by lMoutou et al.l (120041 ) and lGillon et al.l ( l2007l ). 
referred to throughout as 20040517 and 20050420. These light curves contain 274 and 270 
data points, respectively, with samplings of ~60 s and 82 s per exposure (including read- 
out time). The 20050420 dataset includes baseline before and after the eclipse, while the 
20040517 light curve ends right before egress. The average photometric precision of the 
20050420 data is 2.0 mmag per 120 s bin, while the 20040517 data, with more stable pho- 
tometry, has an average dispersion of 0.9 mmag per 120 s time bin. 

Although the original OGLE survey data is much less precise than the other light curves, 
we include it in our analysis in order to test whether its parameters are consistent with the 
more recent data. The OGLE discovery curve, denoted 2002021 9, was observed in the I- 
band and combines a year's worth of full and partial transits (jUdalski et al.l 120031 ). To 
produce the light cu rve shown in Figu re HI we phase-folded the data using the orbital period 
P = 1.689868 days (IGillon et al.ll2007l ). and then binned the data over five minutes intervals 
to decrease the point-to-point scatter before fitting. Due to the faintness of the target and 
the smaller telescope diameter, we note two caveats for the 20020219 transit data: (a) the 
photometric error is much larger than any of the other nine light curves, so this transit has 
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little effect on the parameters that result from a noise-weighted fit (described in the next 
section), and (b) our fitted mi dtime is 6.5 minutes earlier than the value reported (without 
errors) in lUdalski et al.l ( 120031 ). though this discrepancy is also comparable to our revised 
errors, which are about 6 minutes after accounting for correlated noise. 



3. Model and light curve fits 



We fit all ten available light curves jointly, to minimize potential systematics in the 
individual light curves and to derive average parameters of the OGLE-TR-132b transit with 
high precision. We also fit each curve (except the low-precision curve 20020219) individually 
to search for parameter variations over time. 

The light curves we r e fit i n the same manner described in kdams et all J2010al lb[). We 
used the lMandel fc Agoll ( 120021 ) algor ithm to generate analy tical models, and the basic white- 
noise model- fitting code described in lCarter fc WinnI (120091 ) . not inclu ding the wave l et ana l- 
ysis. We used the same limb darkening quadratic law in Equation 1 of Adams et al. ( 2010 al). 
with the initial limb darkening coefficients set to the values computed by I Clare tl (l2000i . 
20041 ) for each filter. The OGLE survey data were observed in /-band, but their precision 
is insufficient to constrain the limb darkening, so we have assumed that the shape of this 
transit follow the same limb darkening law as the Sloan i' transits. The initial values of 
the limb-darken i ng co efficients Ui and U2 for each filter were obtained using th e jktld routine 



bv ISouthworthI (120081 ) . assuming the parameters of the host star derived by iBouchy et al. 
(]2004), = 6411 K, log^ = 4.86, [M/H] = 0.0 and V^ir r n = 2 km s"^ (note that al- 



though these values differ slightly from those of iGillon et al.l ( 120071 ) the slight difference in 
initial limb darkening parameters has negligible impact on the final fitted model parame- 
ters). Throughout we fixed the quadratic term of the limb darkening law, U2] the linear 
term, ui, was also fixed in the individual light curve fits. The orbital period of the planet 
was fixed to P = 1.689868 days, after testing confirmed this parameter has little effect on 
the model, and we assumed throughout that M^, = 1.26 ± 0.03 M(d, -R^, = 1.34 ± 0.08 Rq, 
and Mp = 1.14±0.12 Mj, based on the analysis by lGillon et al.l (120071 ). The free parameters 



were the linear limb darkening coefficient, ui, the planet-star radius ratio, k, the inclination, 
i, the semimajor axis in stellar radii, a/R^, the out-of-transit baseline flux. Foot-, and the 
transit midtime, Tc- 

To determine the best fit value and error of each parame t er, we used a Markov Chain 
Monte Carlo (MCMC) method as described in lCarter h WinnI ( 120091 ). We first fit each light 
curve independently and used the resultant reduced of each fit as a weighting factor for 
each light curve. We then fit all ten light curves simultaneously, assuming common values of 
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k, i, a/R^, and Ui^x (where x indicates the filter), and independent values of Fqot and Tc 
for each light curve. We combined three MCMC chains of 950,000 links each, discarding the 
first 50,000 links of each chain to avoid biasing the solution towards the input values, and 
produced the distributions shown in Figure IH In order to confirm that the chains had truly 
converged, we calculated the Gelman-Rubin statistic for each parameter, and found that it 
was never greater than 1.1, indicating acceptable convergence. 

In Table HJ we report the median value of each parameter distribution, where the error 
bars correspond to the 68.3% credible intervals of the distributions. To account for the effect 
of correlated noise, the errors ha ve been incre a sed ac cording to the residual permutation 
method as previously discussed in lAdams et al.l (l2010al Jbl): with this factor, the noise on our 
transits is 1.2 to 1.9 times greater than the Poisson noise. Table H] also reports two derived 
quantities, the value of the impact parameter, b, and the total duration of the transit, T14, 
defined as the time from the start of ingress to the end of egress. 

To explore potential parameters variations over time, we ran additional MCMC fits 
on the individual transit light curves. We excluded the OGLE survey composite transit, 
20020219, because the photometric precision is much worse than the rest of the data, and 
the model for this light curve is not well constrained without fixing several parameters. 
For the nine light curves we examined individually, the results for four parameters, k, i, 
Ti4 and a/R^f, are summarized in Figure |3l We find no evidence of variation of any of 
those parameters, within the error bars, during the five year time span covered by those 
observations. 

The planetary radius of a transiting planet can be difficult to determine in a crowded 
field, since even very faint comparison stars may dilute the depth of the measured tr a nsit i f 
not carefully accounted for. The two literature light curves published in ICillon et al.l (120071 ) 
were both created using image subtraction, and that work was careful to note that the 
choice of background normalization factor can greatly affect the depth obtained. All seven 
of the new transits were observed during excellent seeing, which allowed us to use aperture 
photometry despite the crowded field. (We note that a star that is about 3.5 magnitudes 
fainter is located 1.3", or 36 unbinned pixels, from our target; this is far enough outside of 
our chosen apertures that its effect would be much less than the millimag precision of our 
light curves.) 

Despite this caveat, we note th at the radius we obtain, Rp = 1.23 ± 0.07 Rj, is very 



similar to t hat of iGillon et al. ( 120071 ). Rp = 1.18±0.07 Rj. It also agrees well with the radius 
derived by ISouthworthI (120081 ) . Rp = 1.25 ± 0.16 Rj, which was deri ved from a reanalysis 
of the higher-preci s ion lig ht curve, 20040517, from iGillon et al.l ( 120071 ). Note that the error 
in the ISouthworthI ( 120081 ) reanalysis, which accounts for correlated noise, is larger than the 
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Gillon et al.l (120071 ) error; this is because the transit is missing the end of egress, which makes 



it particularly vulnerable to correlated noise. With the addition of seven full light curves, the 
planetary radius error can be determined much more precisely, and is currently dominated 
by the uncertainties in the radius of the star, = 1.34 ±0.08 Rq; without the stellar error, 
the radius error would have been 0.01 Rj, seven times smaller. 



4. Timing 

We now turn to an analysis of the timing of each light curve. Each transit midtime, Tc, 
derived in the previous section is listed in Table HI Times for the seven new transits have 
been converted to Barycentric J ulian Days using the B arycentric Dynamical Time standard 



(BJDtdb), as recommended by Eastman et al.l (120101 ). In addition, we have converted the 



three literature datasets to the same BJDtdb time standard by adding the UTC-TT cor- 
rection, after confirming with the authors that none of those data had previously accounted 
for it (M. Gillon and A. Udalski, private communication). 

During the course of this analysis, we noted that the original Tc published for the OGLE 
survey dataset (20020219) does not agree with the more recent transits, prompting us to refit 
the light curve. Our joint fit of all ten transits, which provides high-precision constraints on 
the light curve shape of the OGLE composite curve, gives a midtime of Tc = 2452324. 69684± 
0.00284 BJDtdb, 331 seconds earlie r than the published time (2452324.70067 BJDutc, with 



errors assumed to be about 0.001, lUdalski et al.ll2003l ). After correcting for the UTC-TT 



offset (64.184 s) that timing difference increases to 395 seconds. Our newly calculated transit 
midtime agrees well with the midtimes of the other nine transits, though with quite large 
error bars: the formal error is 245 s, but after accounting for the effects of strong correlated 
noise, we adopt an error of 343 s, very similar in magnitude to the timing shift, although 
the ultimate cause remains unclear. We use our recalculated midtime value for 20020219 in 
the remainder of our analysis, while noting that its large error bars mean it has little effect 
on the weighted fits. 

Assuming a constant period, we derived a new transit ephemeris equation from a 
weighted fit to all ten transits, with the form: 

Tc{N) = 2453142.59199(53) [5 JDrau] + 1.68986531(64)Ar, (1) 

where Tc is the predicted midtime of transit and N is the number of orbital periods since 
the reference midtime. The values in parentheses are the errors in the last two digits. This 
fit, represented in Figure H] by the solid line, has a reduced = 1-2, indicating a good fit; 
the errors reported in the equation have been rescaled slightly, by a factor of v^l.2 = 1.1, to 
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account for the difference from = 1.0. Figure S] shows the timing residuals, or observed 
minus calculated time [O — C), for each transit with respect to the new ephemeris equation 
above. 

The new orbital pe riod of OGLE-TR- 132b is five times more precise than the most 
recent value, reported by lGillon et al.l (120071 ). Such improvement is due to the longer transit 
follow-up time baseline available for this system. Nonetheless, the new orbital period differs 
by less than 0.2 seconds from the one derived by iGillon et al.l ( 120071 ). well within the 1 a 
errors. Additionally, only one individual transit deviates by more than 2 a from the fit 
(20090511, -108 ± 49 s, or 2.2a), with most of the rest within 1 a of the best fit line. There 
is thus no reason to conclude that there are any significant timing variations in this system. 

Using the lack of significant timing variations, we can attempt to constrain the rate of 
orbital decay, and provide upper limits on the masses of potential companion planets. 



4.1. Limits on orbital decay 



Recent calculations by iLevrard et al.l ( 120091 ) predict a lifetime for OGLE-TR-132b of 
only 52 Myr, assuming the tidal dissipation factor for the star of = 10^. (For = 
10^ — 10^°, the lifetime would be 5.2 Myr to 520 Gyr.) Recent observations of a different hot 
Jupiter, 0GLE-TR-113b, hint that its orbital period may be decreasing by Pnsb = — 60± 15 

which, if confirmed, leads to an estimate of = 1.6 x 10^ 
calculations of tidal dissipation. If we assume the same value 



(Adams et al. 


2010b) 


Levrard et al 


( 


2OO9I) 



ms yr 
using the 

of the for OGLE-TR-132, the remaining lifetime of the planet scales to only 0.8 Myr. 

To estimate the predicted rate of decay per year, we define the lifetime of the planet 
as the time that takes for it to reach the Roche lobe of the star. We estimated the radius 



of the Roche lobe using the analytical approximation derived by ( iPaczyhskil Il97ll ) . = 
0.462 a {Mp/{Mp + M^)y/^ = 0.014 AU for OGLE-TR-132b, assuming a circular orbit. The 
corresponding total decrease in period is 1.16 days, so assuming a constant rate of decrease, 
the orbital decay rate is 2 ms yr~^ if = 10^, or a much faster 140 ms yr~^ if = 1-6 x 10^. 
This is, perhaps coincident ally, close to the best fit value of a decreasing-period model if the 
period of OGLE-TR-132b is decaying (Pi32fe = —125 ± 80 ms yr~^, with reduced-x^ = 0.8), 
but with a significance of only 1.5 a this value is entirely consistent with zero period change. 
Since the constant-period ephemeris in Equation [T] fits the data so well, there is no reason 
to conclude that this system is undergoing orbital decay, within the limits of our precision 
and length of observations. 



Additional observations of OGLE-TR-132b over the next decade will be important to 
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refine the limits on orbital decay. 



4.2. Limits on perturber mass 

The lack of significant TTVs for OGLE-TR-132b over the seven year time span of the 
compiled observations indicate the absence of other massive, nearby companions in this sys- 
tem. However, smaller planets can be still hiding within the 50-100 s dispersion of our TTVs. 
In this section we report upper limits to the mass and orbital separation of hypothetical per- 
turbing planets, which would cause TTV signals with amplitudes smaller tha n our detection 



l imit. Those upper limits were computed using an implementation of the ISteffen &: Agol 



( I2OO5I ) algorithm, as implemented by D. Fabrycky (2009, private communication). 



Figure [5] shows the maximum allowed mass for a potential perturbing planet over a range 
of orbital periods. These calculations were performed using only the seven new transits and 
the best literature transit (20040517), due to the higher errors on the other two literature 
points. For simplicity, we assumed that the perturber is coplanar with OGLE-TR-132b, and 
explore several orbital eccentricity values. The most conservative constraints are place by 
objects with the lowest eccentricities, e = and e = 0.05. Objects with higher eccentricities 
(represented here by the e = 0.5 case), might be expected if the ob iect were scattered by 



Kozai migration (e.g.. iFabrycky fc Tremaind 120071 : iTriaud et al.ll2010l ). but near OGLE-TR- 



132b such objects tend to be much more unstable, and to produce larger hypothetical TTV 
signals (and correspondingly lower mass limits). A gray shaded region in Figure [5] shows 
the approximate unstable region for an Earth-mass object on a circular orbit, following 



Barnes fc GreenbergI (120061 ). We examined a range of orbital period ratios from 1/4 and 4 
times that of OGLE-TR-132b (or 0.4 to 7.4 days), divided into 51 bins on a logarithmic 
scale. For each period, we increased the mass of the perturber from an initial value of 0.0001 
Mj, increasing the mass by a factor of 1.1 with each step until we find the mass that would 
produce a TTV equivalent to the 3a confidence level of our results, which we report as the 
upper limit on the perturber's mass. 

The mass constraints placed are strongest at orbital distances near low-order mean mo- 
tion resonances with OGLE-TR-132b. In particular, we can discard the presence of planets 
as small as 5-10 M0 in the external 2:1 and 3:2 and internal 1:2 and 2:3 resonances in cir- 
cular or nearly circular orbits. More eccentric objects have even smaller limits mass limits 
outside of the 1:2 and 2:1 resonances, and are unstable on closer orbits. No significant con- 
straints can be placed on near-circular objects at orbital distances outside the 2:1 resonance, 
at which period perturbers of any mass (including stellar) are essentially unconstrained by 
the timing. Such objects are better constrained by other methods, such as radial velocity. 
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direct detection, and dynamical stability arguments. 

Finally, neither the precision of our light curves in Figure [T] nor the precision of our 
timing residuals in Figure H] are sufficient to constrain the presence of moons around OGLE- 
TR-132b. The best achieved precision on any light curve was about 1 mmag on the photo- 
metric depth, 37 s on the transit midtime, and 213 s on the transit duration, whereas an 
Earth-sized planet located within the Hill radius of the planet would produce transit signals 
of 0.05 mmag, a TTV smaller than a few seconds, an d a transit duration variation, or TDV, 
of a few seconds (lKippingll2009l : iKipping et al.ll2009l ). 



5. Summary and Conclusions 



We have measured seven new transits of the hot Jupiter OGLE-TR-132b, which triples 
the number of transits available for this planet to date. The new transits have timing 
precisions between 47 and 87 seconds and photometric precision from 1-2 mmag in two 
minute bins. The updated observational dataset of OGLE-TR-132b spans seven years (2002- 
2009) and a total of ten transits after combining the new transits presented in this paper 
with three previously published transit epochs. 

We find no evidence of transit timing variations with amplitudes larger than —108 ±49 s, 
which allows us to place significant constraints to the presence of additional close-in planets 
in the system. In particular, near both the internal and external 2:1 and 3:2 mean motion 
resonances, objects as small as 5-10 Earth masses would have been detectable if they existed. 

We also tested for orbital period decay, with the best fit rate of period decrease, P1326 = 
— 125 ± 80 ms yr~^, consistent with no change. The predicted orbital period decay rate for 
this planet is 20 ms yr~^ if = 10^, and 140 ms yr~^ if Q^, = 1.6 x 10"^ as suggested for the 
star OGLE-TR-113 (lAdams et al.ll2010bl ). so additional transits with better timing precision 
(< 30 s) or a significantly longer time baseline are needed to search for slower rates of period 
change. 

Finally, we have measured a revised radius val ue for the planet, R p = 1.23 ± 0.07 Rj, 



whi ch is fully co n sisten t wit h the radius derived by ISouthworthI (120081 ) from a reanalysis of 
the iGillon et al.l (120071 ) and iMoutou et al.l (120041 ) data. We used simple aperture photom- 
etry to extract all transit light curves, so the depth of our transits should not be affected 
by the potential normalization pro blems found when using image differencing techniques, 
as reported for this same planet by iGillon et al.l (120071 ). Our planet ra dius provides an in- 
dependent confirmation, within errors, of the transit depth obtained by IGillon et al.l (120071 ) 
using image deconvolution techniques. 
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Fig. 1. — Ten transits of OGLE-TR-132b. Normalized flux for all available light curves are 
plotted vs. orbital phase, with each transit shifted vertically by 0.015 for clarity. Though 
the full data was used for fits, the data displayed has been binned to 2 minutes (except for 
the OGLE composite c urve 20020219 , whic h wa s fit with 5 minute bins). The top three 
transits are taken from lUdalski et al.l (|2003[ ) and iGillon et al.l ( 120071 ). while the rest of the 
transits are new. The joint model fit (solid lines) were calculated using the parameter values 
in Table IH the stated standard deviations are the residuals of each transit from the joint 
model fit. 
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Fig. 2. — Parameter distributions for joint fit to ten transits of OGLE-TR-132b. Smoothed 
histogram of normahzed parameter distributions, from which the parameters in Table H] are 
derived. The sohd hne is the median value (which is very close to the mean value in all 
cases). The dashed lines show the 68.3% credible interval. These values were calculated for 
2,850,000 links. 
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Fig. 3. — Parameter variation of nine individual transits of OGLE-TR-132b, based on indi- 
vidual MCMC fits (Table E]), for (clockwise, from top-left): k, i, a/R^, and T14. The OGLE 
composite curve 20020219 has been omitted since it was not fit independently. The values 
derived from the joint MCMC fit to all transits (Table H]) are plotted as solid black lines with 
±1(T errors. The dashed red lines indicate the best sloped fit with ±la errors. All fits are 
consistent with no parameter variation during the time period examined. 
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Fig. 4. — Observed minus calculated midtimes for OGLE-TR-132b. Timing residuals using 
the new ephemeris (Equation [1]), calculated using all ten transit epochs. The solid line 
represents zero deviation from expected time of transit, while the dashed lines represent the 
la and 3a errors on the calculated orbital period, indicating the slopes that result for a 
mis-determined period. 
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Fig. 5. — Upper mass limit on additional planets with initial eccentricities = 0.0 (red), 
Cc = 0.05 (black) and Cc = 0.5 (gray). The central point shows 0GLE-TR-113b. Near the 1:2 
and 2:3 mean-motion resonances, objects as small as 5-10 would have been detectable. 
The shad ed grey region shows the ins tability region for a 1 companion in a circular orbit, 
following iBarnes fc GreenbergI (120061 ) . 
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Table 1. Observational parameters for seven transits of OGLE-TR-132b 



Transit 


Frames'' 


Exp. Time 


Binning 


Readout 


Airmass 


Time range observed 


Atm. stability 


Seeing 


(UT) 




(sec) 




(sec) 










20080419 


227 (0) 


60 


1x1 


5 


1.2-1.4 


4.5 hrs around midtraiisit 


Stable 


0."4-0."5 


20080424 


160 (0) 


60, 120 


1x1 


5 


1.2-1.8 


3 hrs (25 min before ingress) 


Variable 


0."6-0."7 


20080511 


136 (5) 


120 


1x1 


5 


1.2-1.5 


5 hrs around midtransit 


Stable 


0."5-0."7 


20090207 


1190 (10)"= 


10, 15 


2x2 


d 


1.2-2.1 


5 hrs around midtransit 


Variable 


0."4-0."7 


20090311 


2322 (0) 


10, 15 


2x2 


d 


1.2-1.6 


8 hrs (4 hrs before ingress) 


Very stable 


0."6 


20090424 


476 (0) 


20-75 


1x1 


d 


1.2-1.6 


5 hrs (30 min before ingress) 


Variable 


0."3-0."4 


20090511 


440 (0) 


35, 40 


2x2 


d 


1.2-1.8 


5 hrs (35 min before ingress) 


Stable 


0."5 



^Number of frames used (additional frames that were discarded). 
''Discarded due to tracking problem. 
'^Discarded due to low counts. 
Readout is a few miliseconds in frame transfer mode. 



Table 2. Photometry parameters for OGLE-TR-132b 



Transit 


Comp. Stars 


Aperture*^ 


Sky radius, width 


Slope removed 


Scatter'' 


Noise factor'^ 


(UT) 




(pixels) 


(pixels) 




(mmag) 




20080419 


7 


13 


30, 20 


None 


1.5 


1.3 


20080424 


4 


19 


40, 30 


+0.013% d 


2.2 


1.9 


20080511 


3 


26 


50, 20 


-1.2% " 


1.3 


1.4 


20090207 


2 


15 


30, 10 


+0.1% ° 


3.9 


1.6 


20090311 


2 


12 


30, 10 


-0.3% ° 


2.9 


1.2 


20090424 


3 


15 


30, 10 


-0.052% f 


2.0 


1.2 


20090511 


8 


9 


20, 10 


+0.12% 


1.8 


1.2 



''Radius around star 

''Prame-by-frame scatter on out-of-transit flux (not corrected for exposure time). 

•^Factor by which scatter exceeds Poisson noise estimate. 

'^Trend removed against telescope azimuth, in units of flux deg~' . 

"Trend removed against object airmass, in units of flux X~'. 

'Trend removed against seeing in pixels, in units of flux pixels"'. 
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Table 3. Flux values for seven transits of OGLE-TR-132b 



Mid-exposure (UTC)'' Mid-exposure (BJD) Normalized stellar flux RMS error 

2454575.506679 2454575.510331 0.9996423 0.00122 

2454575.507433 2454575.511084 1.001028 0.00122 

2454575.508183 2454575.511834 0.9992174 0.00122 

2454575.508936 2454575.512588 1.000828 0.00122 

2454575.509689 2454575.513341 1.000836 0.00122 



^This table is available in its entirety in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content. 



Table 4. Transit parameters for OGLE-TR-132b (jointly fit) 



Parameter Median value'' 



Fitted Parameters 



k 




0.0941 ±0.0014 




a/R, 




4.81 ±0.16 




i (dej 


i) 


83.6 ±0.7 








0.24 ± 0.04 








0.19 (fixed) 








0.23 ±0.06 








0.38 (fixed) 




Tc- 


2452324 


0.69684 ± 0.00400 


(343 s) 


Tc- 


2453142 


0.59191 ± 0.00043 


(37 s) 


Tc- 


2453480 


0.56473 ±0.00110 


(94 s) 


Tc- 


2454575 


0.59877 ± 0.00074 


(64 s) 


Tc- 


2454580 


0.66772 ± 0.00100 


(87 s) 


Tc- 


2454597 


0.56611 ±0.00054 


(47 s) 


Tc- 


2454869 


0.63518 ± 0.00055 


(47 s) 


Tc- 


2454901 


0.74172 ±0.00079 


(68 s) 


Tc- 


2454945 


0.67790 ± 0.00075 


(64 s) 


Tc- 


2454962 


0.57567 ± 0.00057 


(49 s) 



Derived Parameters 

b 0.53 ± 0.05 

Ti4 (sec) 9322 ± 111 

Rp (Rj)^ 1.227 ±0.074 

a (AV)^ 0.030 ± 0.002 



^Median value of parameter distribution with errors in- 
creased from the 68.3% credible interval by residual permu- 
tation method (see Section [S]!. 

''Assuming i?« = 1.34 ± 0.08 Rq jCillon et al.ll2007l ') and 
using Rj = 71, 492 km. 
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Table 5. Individual transit parameters for OGLE-TR-132b (independently fit) 



Transit 






Ti4^ 












20040517 


0.0933 ± 0.0016 




9310 ± 213 


1.1 


4.86 ± 0.36 




83.9 ± 1.9 




20050420 


0.0938 ± 0.0036 




9627 ± 654 




4.50 ± 0.78 




82.1 ±4.0 




20080419 


0.0882 ± 0.0043 


4.1 


9041 ± 300 


2.9 


5.46 ± 0.42 


1.9 


87.6 ± 1.9 


1.2 


20080424 


0.1002 ±0.0027 




9703 ± 384 


1.1 


4.25 ± 0.57 


1.4 


80.8 ± 2.6 


1.4 


20080511 


0.0955 ±0.0024 


1.6 


9305 ± 226 


1.3 


4.97 ± 0.35 




84.6 ±2.0 




20090207 


0.0939 ± 0.0016 




8919 ± 223 


1.3 


5.35 ± 0.40 




86.0 ±2.3 




20090311 


0.1014 ±0.0053 


4.6 


9523 ± 245 


1.6 


4.55 ± 0.23 




82.3 ± 1.1 




20090424 


0.0918 ± 0.0029 


1.7 


9286 ± 228 


1.1 


4.76 ± 0.40 




83.3 ±2.1 




20090511 


0.0911 ± 0.0026 


1.6 


9150 ± 214 


1.0 


4.96 ± 0.42 




84.2 ± 2.2 





''Formal individual MCMC fit value and error (scaled upward by factor / in adjacent column). 

''Factor by which the error in the previous column has been increased based on the residual permu- 
tation method; no value is given if the error from residual permutation is smaller than the MCMC fit 
error. 



Table 6. Timing residuals for OGLE-TR-132b 



Transit Number O - C (s) ao-c/{0-C) 



20020219 


-484 


-29 ± 343 


-0.1 


20040517 





-7 ±37 


-0.2 


20050420 


200 


-28 ± 94 


-0.3 


20080419 


848 


87 ±64 


1.4 


20080424 


851 


31 ±87 


0.4 


20080511 


861 


8 ±47 


0.2 


20090207 


1022 


73 ±47 


1.6 


20090311 


1041 


-5 ±68 


-0.1 


20090424 


1067 


-32 ± 64 


-0.5 


20090511 


1077 


-108 ±49 


-2.2 



